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Spectral Test of the Goodness of Fit of ARIMA Model 


WSTĘP 


Widmo mocy procesów stochastycznych stanowi ważne narzędzie 
w analizie zjawisk ekonomicznych. Pozwala ono bowiem określić wpływ 
różnych częstotliwości na wariancję procesu oraz wyróżnić istotne waha- 
nia periodyczne. Badania, których rezultaty zawarto w pracy ([4]) po-- 
twierdzają hipotezę o celowości wykorzystania funkcji spektralnej do 
oceny wariancji procesu z istotną autokorelacją. 

Celem niniejszej pracy jest zastosowanie funkcji widmowej do oceny 
zgodności modelu ARIMA !, z danymi empirycznymi. 

Pierwsza część rozważań poświęcona jest zagadnieniom teoretycznym 
związanym z procesami ARIMA i testem równoważności widma mocy. 
Rezultaty identyfikacji modelu opisującego kształtowanie się cen giełdo- 
wych złota w latach 1977—1979 oraz ocenę zgodności tego modelu przy 
pomocy proponowanej przez autora procedury zawarto w końcowej części 


pracy. 


LINIOWE MODELE TYPU ARIMA I ICH IDENTYFIKACJA 


Liniowe modele typu ARIMA wykorzystuje się do opisu zarówno sta- 
cjonarnych, jak i niestacjonarnych procesów losowych. Niniejsze rozwa- 
żania ograniczono jedynie do tych procesów, które okazały się istotne 


1 Ang. Autoregressive Integrated Moving Average. 
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z punktu widzenia budowy modelu przedstawiajacego ksztattowanie sie 
cen giełdowych złota w latach 1977—1979. Dokładną prezentację szero- 
kiej klasy procesów ARIMA zawarto m.in. w pracach ([1]) i ({3}). 
Niech B oznacza operator cofnięcia 2, zaś V operator różnicowania 3 o pod- 
stawie 1, czyli: 


Bz, = 2-1 (1) 
Bmz, = Z—-m (2) 
Vz = Z:— 241 = (1—B) z; (3) 


gdzie (z+) oznacza ciąg losowy. 


Przy takich założeniach ogólną postać procesu typu ARIMA (p,d,q) 
można przedstawić w następujący sposób: 


F(B) Vaz, = G(B) a (4) 
gdzie 
{a} — proces czysto losowy 4 
{4} — badany ciąg losowy 
F(B) = 1-1,B—1,B?— ... —1,BP 


Viz = (1—B)d z, = we 


W przypadku, gdy d= 0 proces (4) zaliczany jest do klasy procesów 
ARMA (p, q). Jeżeli d=0 i q=0 wzór (4) odpowiada procesowi auto- 
regresyjnemu rzędu p(AR(p)). Dla d=0 i p=0 mamy do czynienia 
z procesem średniej ruchomej rzędu q (MA(q)). Operacja różnicowania 
Ndz, ma na celu doprowadzenie ciągu (z+) do stacjonarności.5 

Jednym z najistotniejszych zagadnień związanych z procesami auto- 
regresyjnymi jest identyfikacja modeli opisujących te procesy. Wydaje 
się, że z praktycznego punktu widzenia najbardziej trafne uwagi zawarto 
w pracy (3, s. 173—207). Ogólnie mówiąc zagadnienie identyfikacji można 
sprowadzić do następujących etapów: 

1) takiego różnicowania ciągu {z}, które dałoby stacjonarny 
ciąg (wz). | 

2 Ang. backward shift operator. 


3 Ang. backward difference operator. 
4 Proces jest czysto losowy, gdy: 


E(a+) = o 
: o s#0 
i Eat, at+s) = ee <= 0 


5 Autor ma na myśli stacjonarność w szerszym sensie. Ciekawe uwagi odnośnie 
poszczególnych rodzajów stacjonarności zawarto m.in. w pracy (11). 
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2) identyfikacji modelu typu ARMA dla {w;}. 

Podstawę procedury  identyfikacyjnej zaproponowanej przez 
G.E.P. Boxa i G. M. Jenkinsa stanowi zachowanie się funkcji autokore- 
lacji (ox) i autokorelacji cząstkowej (lyk) procesu € {w,}. Jeżeli ciąg {w} 
uzyskany z {z} przez różnicowanie (w, = Vaz) charakteryzuje się dość 
szybko wygasającą funkcją autokorelacji, to istnieją przesłanki do stwięr- 
dzenia, że (w,) jest stacjonarny. 

Można więc w takim przypadku przystąpić do identyfikacji procesu 
ARMA (p,q). Autorzy pracy (3, s. 175) stwierdzają, że dla procesu mie- 
szanego zawierającego AR(p) i MA(q) w zakresie funkcji autokorelacji 
istnieje możliwość zaobserwowania nakładającej się gasnącej fali sinuso- 
idalnej i wykładniczej po pierwszych (q-p) rzędach. W zachowaniu się 
funkcji autokorelacji cząstkowej podobną prawidłowość dostrzega się po 
(p-q) rzędach. Przedstawione w końcowej części pracy badania empirycz- 
ne wymagają krótkiej charakterystyki procesów ARIMA o następujących 
parametrach: a) (2, 2, 0), b) (0, 1, 0), c) (0, 1, 1). 

Ad. a) Proces ARIMA (2, 2, 0) wyraża się zależnością postaci: 

(1—1,B—1,B?) V/2z, = a, (5) 
gdzie 
{a} jest procesem czysto losowym. 


Załóżmy, że r; (1i=1,2,..,k) oznacza ocenę teoretycznej autokorelacji 
rzędu (i) ciągu V?2z: (t =3,....n). Oceny parametrów 1, i lą można wy- 
znaczyć na podstawie następujących równań: 


~ ry (1-15) 
> ZE ee 6 
h 1 r? > ( ) 
A Ia —r? 
= 7 
le 1 —r? ( ) 


Macierz wariancji i kowariancji estymatorów l; i lą szacujemy wykorzy- 
stując wzór postaci: 
ee 1-5 —i (+h) 
V (le) = 07? | A ź a2 (8) 
—h(l+h) 1-5 
Podstawę estymacji wariancji procesu czysto losowego {a,} stanowi na- 
stępująca zależność: 


83 = c, (—rlh —Tal>) (9) 
gdzie 
5,2 — ocena wariancji (a), 
co — ocena wariancji {w+}. 


« Dokładne omówienie tych funkcji można znaleźć w pracach: (1) i (3). 
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W pracy (3, s. 176) podkresla sie, ze proces jest typu ARIMA (2, 2, 0), 
jeśli wykres funkcji autokorelacji ciągu {wi} (We = Zt 24-1 t 2-2; 
t= 3, ...,n) można uzyskać z nakładającej się krzywej wykładniczej i tłu- 
mionej fali sinusoidalnej. Funkcja autokorelacji cząstkowej takiego pro- 
cesu powinna przyjmować tylko dwie pierwsze różne od zera wartości. 

Ad. b) Jeżeli okaże się, że funkcja autokorelacji ciągu {wi} (wi = 
= W4 Z1 ; t=2,...,n) charakteryzuje się wykresem zbliżonym do wy- 
gasającej funkcji wykładniczej, zaś wszystkie wartości funkcji autokore- 
lacji cząstkowej nie różnią się istotnie od zera; to będą podstawy do opi- 
su procesu {z,} przy pomocy modelu: 


Vz=a (t=2,3,..,n) (10) 
gdzie 
{ay} jest procesem czysto losowym 


Wariancje (a+) można oszacować wyznaczając ocenę tego parametru 
dla (w). 

Ad. c) W przypadku, gdy funkcja autokorelacji procesu {wi} (wą = 
= Vz) przyjmuje wartość różną od zera jedynie dla pierwszego rzędu, 
zaś funkcja autokorelacji cząstkowej charakteryzuje się przebiegiem zbli- 
żonym do wykładniczego, istnieje możliwość wysunięcia hipotezy, że pro- 
ces jest typu ARIMA (0,1, 1). Proces ten wyraża się związkiem postaci: 


Vz =a h'a  (£=2,38,.., n) (11) 
Parametr h, szacujemy rozwiązując następujące równanie kwadratowe: 
—h 
= R 12 


Sredni błąd szacunku h; uzyskuje się z zależności: 
FK) =n"' (1—h$), (13) 
zaś wariancję procesu (a+) wyraża się wzorem: 
52 =c,:(1-i-hq). (14) 


Na zakończenie rozważań odnośnie procesów ARIMA warto wspomnieć 
o zagadnieniu weryfikacji hipotezy o zgodności danego modelu z danymi 
empirycznymi. W większości prac z omawianego zakresu szczególną rolę 
przywiązuje się do badania własności reszt odpowiadających określonym 
modelom. Do weryfikacji postawionej powyżej hipotezy wykorzystuje się 
statystykę: 
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Q=n)rffa), (15) 
k=1 
gdzie i 
r,(a) — ocena wartości funkcji autokorelacji reszt modelu ARIMA, 


która przy założeniu prawdziwości hipotezy zerowej ma rozkład chi — 
kwadrat z(m—p—q) stopniami swobody. Jeżeli wartość statystyki Q 
okaże się wyższa od wartości krytycznej Xj, to hipotezę o zgodności 
modelu z danymi empirycznymi należy odrzucić. 

W niniejszej pracy zaproponowana zostanie procedura pozwalająca zwe- 
ryfikować hipotezę zerową na podstawie porównania widm mocy proce- 
cu teoretycznego i rzeczywistego. 


TEST RÓWNOWAŻNOŚCI WIDMA MOCY W OCENIE ZGODNOŚCI MODELU 
AUTOREGRESYJNEGO Z DANYMI EMPIRYCZNYMI 


W pracach z zakresu analizy spektralnej zwraea się uwagę na 
wszechstronne możliwości wykorzystania funkcji widma mocy do bada- 
nia własności stacjonarnych ciągów losowych. Wśród licznych zastosowań 
powyższej funkcji warto podkreślić jej następujące zalety: 

a) informuje o wpływie poszczególnych częstotliwości na wariancję 
procesu, 

b) pozwala wyznaczyć w danym ciągu losowym istotne wahania pe- 
riodyczne, 

c) może stanowić podstawę weryfikacji hipotezy, że dany proces jest 
czysto losowy. f 

W pracy (11) widmo mocy procesu zostało wykorzystane do wyodręb- 
nienia jednorodnych podzbiorów z punktu widzenia określonej cechy. 
Autorzy prac (3, s. 294—295) i (8, s. 334—359) zwracają uwagę na możli- 
wość zastosowania omawianej funkcji do oceny zgodności modelu ekono- 
metrycznego z danymi empirycznymi. Przy założeniu, że ciąg (a+) repre- 
zentuje reszty modelu autoregresyjnego, stopień dopasowania tego mode- 
lu do procesu rzeczywistego można określić na podstawie wykresu 


funkcji C(fj) postaci: 
j 
Sri) 


C(fj) = aS (16) 
gdzie 7 
I(fi) = 2n7! | (a, cos 2nfit)? + Ya, sin ati 


fi=in; (i=0,1,... [n/2]8) 
ôa? oznacza ocenę wariancji procesu {a}. 


12 — Annales 
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Jak wykazują G.E.P. Box i G.M. Jenkins (3, s. 295) analizowany 
model można uważać za dobrze dopasowany, jeżeli wartości funkcji C(fj) 
oscylują% wokół prostej przechodzącej przez punkty o współrzędnych 
(0, 0) i (0,5, 1,0). | i 

W niniejszym opracowaniu proponuje się metodę, która jak się wy- 
daje powinna opisywać postawiony problem w sposób bardziej jedno- 
znaczny. Załóżmy, że rozważamy proces losowy typu ARIMA (p, d, q) 
postaci: 


(1—1,B—1,B?— ... —1,BP) w, = (1—h;B— ... —hgB9) at (17) 
gdzie i 
Wt = Wt Ww 
w; = Viz 
1 n=d ; 
l; (i= 1, 2, ..., p) — oceny parametrów procesu autoregresyjnego, : 
h; (i = 1, 2,..,q) — oceny parametrów procesu średniej ruchomej, 
© {a} — proces czysto losowy. 


Jak łatwo zauważyć proces wyrażający się zależnością (17) można przed- 
_ stawić w postaci: 


wę = (4B+1,B?+ ... +1,BP) w; +(1—h;yB—hgB?— ... —hgB9) a, (18) 


Dla realizacji niestacjonarnego procesu (z+) utwórzmy następujące ciągi 
losowe: 


-e ks 
w =W 


wr = lW, ;-+...-Hlow, _ -+a,—hya,_,—...—h,a,_, (19) 


P 
gdzie 
wę = Viz- W 
ay; (t= 1,2, ...,n-d) — ciąg liczb losowych uzyskanych przy wykorzy- 
staniu generatora liczb losowych o rozkładzie 
N(0, 6,2). 


Postawmy hipotezę, ze model ARIMA (p, d, q) jest dobrze dopasowa- 
ny do danych empirycznych. Wydaje się, że przy założeniu prawdziwości 


7 Wykres I(fi) w zależności od fi nazywany jest periodogramem, 

8 [m/2] oznacza część całkowitą z liczby n/2. 

s Słowo „oscylują” stwarza możliwość subiektywnego spojrzenia na problem 
zgodności modelu. 


D 
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powyższej hipotezy ciągi: (w;T) i (w;*) powinny reprezentować ten sam 
proces. Widma mocy ! tych ciągów nie mogą więc różnić się od siebie 
w sposób statystycznie istotny. Załóżmy, że Ho ma postać: 


Pr(f) = Pelf) 


gdzie 


Pr(f) — widmo mocy procesu (w;T), 
Pz(f) — widmo mocy (w+F). 


Weryfikację hipotezy Ho można przeprowadzić na podstawie testu rów- 
noważności widma mocy zaproponowanego przez J. S. Bendata 
i A. G. Piersola w pracy (2, s. 313).-Autorzy testu wykazali, że przy za- 
łożeniu prawdziwości H, statystyka Dgr? wyrażona wzorem: 


2 PAŃ W Pre |? 
2 2, A eee be ame 
Diz (s; i = > [is PEC) (20) 
gdzie 
Pz(ij) — oceny widma mocy odpowiednio procesu {w,T} i {wi}, 


(m+ 1) — liczba punktów estymacji funkcji spektralnej, 
Sk, Sr — liczba stopni swobody estymatorów Py(f) i Pr(f), 
ma rozkład chi-kwadrat z (m+ 1) stopniami swobody. 


Jeżeli okaże się, że wyznaczona wartość D2, spełnia nierówność 
Dł, >X? (X2 jest wartością krytyczną odczytaną z tablicy rozkładu 
chi-kwadrat na poziomie istotności a i dla (m+ 1) stopni swobody), hipo- 
tezę o równoważności badanych widm należy odrzucić. Istnieją więc 
w takim przypadku podstawy do stwierdzenia, że model nie jest zgodny 
z danymi empirycznymi. 


OCENA ZGODNOŚCI AUTOREGRESYJNEGO MODELU OPISUJĄCEGO 
KSZTAŁTOWANIE SIĘ CEN GIEŁDOWYCH ZŁOTA W LATACH 1977—1979 


Podstawę rozważań tej części pracy stanowi szereg tygodniowych no- 
towań cen złota na giełdzie w Londynie w latach !! 1977—1979. Szereg 
ten przedstawiono na ryc. 1. 


10 Pod pojęciem widma mocy stacjonarnego procesu losowego rozumie się trans- 
formatę Fouriera jego funkcji kowariancji. i 

11 Dane te zaczerpnięto z tygodnika „Życie Gospodarcze” (poszczególne numery 
w latach: 1977, 1978, 1979). ges) 


127 


180 Adam Góral 


Ryc. 1. Ceny złota na giełdzie w Londynie w latach 1977—1979 (w dol. za troy uncję) 
Prices of gold at the Exchange in London in the years 1977—1979 (in dollars for 
a troy ounce) 


Już pobieżna analiza powyższego wykresu wskazuje na niestacjo- 
narność omawianego procesu. Postawiono więc hipotezę, że opis procesu 
generującego ceny giełdowe złota (w dol. za troy uncję) w Londynie 
w latach 1977—1979 może być dokonany na podstawie modelu typu 
ARIMA (p, d, q). W pierwszym etapie identyfikacji przeprowadzono ope- 
rację różnicowania na wyjściowym ciągu {z} (t= 1,2,..., 155). Wartości 
funkcji autokorelacji i autokorelacji cząstkowej rzędu 1—20 oszacowano 
dla następujących ciągów: 


z, = ZZ, 
—w (i=1,...,53 t =i+1,...n) 
gdzie 
wi = F'z, 
Pierwszą z wymienionych funkcji estymowano na podstawie wzoru 
postaci: 
Tk = C/Co ` (21) 


gdzie !? 


n—k 
) WiWiie 
_ ti 


n 


Cy 3 


12 Jest to asymptotycznie nieobciążony estymator funkcji autokowarigncji pro- 
cesu losowego [1, s. 477]. 
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n 
(E ) w? 
t=1 


n — liczba elementów ciągu {w;}. 


Oceny funkcji autokorelacji cząstkowej wyznaczono rozwiązując układ 
równań Yule'a-Walkera, który w zapisie macierzowym można przedsta- 
wić równością: 


Ri=r (22) 
gdzie 
R — macierz ocen wartości funkcji autokorelacji: 
1 Try Ta Ty 
= rj 1 Iz Ik-2 
Ik-1 Tk-2 Tk-1 l 
1 — wektor ocen współczynników autoregresji procesu AR(k): 
Ik, 
i=|* 
lee, 
r — wektor ocen funkcji autokorelacji 
= 
r=| > 
ko 


Wartości 13 lkk dla k=1,2,..., L (w naszym przypadku L=20) stanowią 
oceny funkcji autokorelacji cząstkowej. Rezultaty uzyskane w przedsta- 
wiony w powyżej sposób zawarto w tabelach 1 i 2. 

Uważna analiza wartości r oraz l (i= 0, 1,..,5) w konfrontacji 
z uwagami zawartymi w części 2, pozwoliła wysunąć przypuszczenie, że 
do badanego szeregu można dopasować następujące modele typu ARIMA: 
(2, 2,0), (0,1,0) i (0,1,1). Do szacowania parametrów występujących 
w pierwszym i trzecim modelu wykorzystano wzory: (6), (7) i (12). Śred- 
nie błędy szacunku tych parametrów uzyskano natomiast na podstawie 
zależności: (8) i (13). 


—r2 
T Ti 


—r2 
l-r 


13 Łatwo zauważyć, że np. lu =, l= 
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Tab. 1. Wartości estymatorów autokorelacji (rx) ~ 
Values of the estimators of the autocorrelation (rx) 


Postaé ciągu losowego 


ka ; 
Zt Ve V2z V 32 Vaze VSzę 
1 0,9605> —0,1780 —0,6184 —0,7158  —0,7616 —0,7969 
2 0,9297 0,1084 0,0958 0,2077 0,3048 0,3909 
3 0,8973 0,1218 0,1047 0,0607 0,0066 0,0787 
4 0,8716 —0,0602 —0,0683 —0,0688  —0,0634 —0,0533 
5 0,8476 —0,0702 —0,0032 0,0412 0,0808 0,1189 
6 0,8230 —0,0998 _ —0,0858 —01237  —0,1719 —0,2119 
7 0,8031 0,0789 0,2155 0,2678 0,2967 0,3055 
8 0,7834 —61965  —0,3206 —0,3508 _ —0,3497 —0,3258 
9 0,7611 0,2586 0,2973 0,3022 0,2802 0,2534 
10 0,7313 —0,0491 —0,1456 —0,1485  —0,1423 —0,1387 
11 0,7105 5,1219 0,0137 —0,0181 0,0276 0,0364 
12 0,6772 0,1017 0,0436 0,0423 0,0417 0,0397 
13 0,6494 —0,0029  —0,0615 —0,0738 _ —0,0827 — 0,0861 
14 0,6288 + 0,0468 0,0688 0,0920 0,1032 0,1081 
15 0,6137 —0,0469 _ —0,0995 —0,1087 _ —0,1149 —0,1176 
16 0,6031 0,1125 0,1007 0,1101 0,1159 0,1143 
17 0,5905 0,0176 0,0742 —0,1002 —0,1051 —0,0980 
18 0,5809 0,1239 0,1007 0,0992 0,0901 0,0799 
19 0,5681 —0,0401  —0,0880 —0,0770  —0,0662 —0,0620 
20 0,5606 0,0678 0,0211 0,0214 0,0248 0,0319 


a k — oznacza rząd autokorelacji, 
b wartość ta jest oceną autokorelacji rzędu pierwszego ciągu {zt}. 
Żródło: Obliczenia własne. 


Tab, 2. Wartości estymatorów autokorelacji cząstkowej (lx) 
Values of the estimators of partial autocorrelation (Txx) 


Postać ciągu losowego 


ka 
Zt Vze V2z+ V3z V4zę VSzę 
1 —0,9695> —0,1780 —0,6184 —0,7158 — 0,7616 —0,7969 
2 —0,0943 0,0793  —0,4641 —0,6248  —0,6557  —0,6689 
3 —0,0273 0,1597 —0,1971 —0,4877 —0,5082 —0,4554 
4 0,0673  —0,0231  —0,0521 —0,3586  —03708 _ —0,2915 
5 0,0280 —0,1224 —0,0206 —0,0995 0,0564 0,1857 
6 —0,0159  —01538  —0,2452 —0,3299 —01287  —0,0162 
7 0,0522 0,0793 0,0069 —0,0455 0,1037 0,0755 
8 "0,0107  —0,1218  —0.2866 —0,2468  —0,0327  ~—0,0143 
9 — 0,0426 0,2369 — 0,0635 — 0,1931 — 0,0579 — 0,0136 
10 — 0,1094 0,0230 —0,0770 — 0,1389 —0,0714 —0,0151 
11 0,0891 0,1142  —0,0381 0.0041  —0,0330  —0,0368 
12 —0,1597 0,0518 —0,0484 0,1894 0,0902 0,0222 
13 0,0054 0,0085 — 0,0966 0,2632 0,1395 0,0916 
14 0,1114 —0,0029 —0,0979 0,2595 0,0573 0,0839 
15 0,0563 0,0138 —0,1033 0,2510 —0,0552 — 0,0052 
16 0,0475 0,0959 — 0,1492 0,1582 —0,1171 — 0,0561 
17 0,0155 0,1900 — 0,0959 — 0,0439 — 0,1555 — 0,0931 
18 0,0406 0,1255 0,0401 — 0,1196 — 0,0092 0,0344 
19 — 0,0309 — 0,0203 0,0184 —0,2158 0,0701 0,0517 
20 0,0692 — 0,0232 — 0,0024 — 0,3318 0,0746 0,0087 


a k — oznacza rząd autokorelacji cząstkowej, 
b wartość ta oznacza ocenę autokorelacji cząstkowej rzędu pierwszego ciągu {zı}. 


Źródło: Obliczenia własne, 
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Efektem prowadzonej w ten sposób estymacji były modele postaci: 


a) w, = —0,9-w,_, —0,45-w,_ +a: (23) 
(0,0012) (0,0012) 

b) w, =a (24) 

c) w, =a,—0,19a,_, (25) 


Kolejny etap badania stanowiła weryfikacja hipotezy o zgodności pro- 
ponowanych modeli z danymi empirycznymi. Weryfikacji tej hipotezy 
dokonano dwoma metodami. Pierwsza z nich wykorzystywała wyrażoną 
wzorem (15) statystykę Q, zaś druga statystykę D2, (20). Wyznaczenie 
wartości Dh, wymagało oszacowania funkcji widma mocy procesów: 
{w,T} i {wE}. Ze względu na niezbyt długie szeregi czasowe (155 obser- 
wacji) oraz konieczność wcześniejszego wykorzystania ocen funkcji auto- 
kowariancji zdecydowano się szacować funkcję spektralną metodą stan- 
dardową !4 wykorzystując wzór postaci: 


plf) =2 | co+2 È` 1/2 (1 +c0s (zh/m)) -ca cos xt (26) 
h=1 
gdzie 
p(fx) — ocena funkcji widma mocy dla fy = k/2m (k = 
= 0, 1,..., m), 
Ch; (h = 0, 1, 2,...., m) — ocena funkcji autokowariancji rzędu h. 


Szczegółowe uwagi odnośnie tego estymatora zawarto w większości prac 
z zakresu analizy spektralnej 5. W niniejszych badaniach zgodnie z suge- 
stią K. Grangera i M. Hatanaki (5) wartość punktu odcięcia funkcji auto- 
kowariancji przyjęto na poziomie m = n/5. 

Liczbę stopni swobody estymatorów Px(fx) i Pr(fx) wyznaczono z za- 
leżności: 


Sr = Sg = 2 ° n/m 


Wartości statystyk: Q i Dgr? dla modeli (23), (24) oraz (25) wraz z od- 
powiednimi wartościami krytycznymi 16 zawarto w tabeli 3. 


“u W tym samym celu można było również wykorzystać metodę szybkiej trans- 
formacji Fouriera (FFT). 

15 Zob.: (2), (5), (7), (9), (10). r 

ie Wartości te odczytano z tablicy rozkładu chi-kwadrat dla poziomu istotno- 
ści a = 0,01 oraz odpowiedniej liczby stopni swobody. 
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Tab. 3. Ocena zgodności modelu z danymi empirycznymi 
The estimation of the agreement of the model with empinical data ` 


Postać modelu Wartość Q xXx? Wartość x? 
28,1615 37,1564 189,0380 52,08 
42,1779 38,5822 56,7838 52,0 
33,8498 38,5822 47,3864 52,0 


a Wartość te uzyskano poprzez interpolacje liniową wartości x2 dla .30 i 40 
stopni swobody (a = 0,01). 


Źródło: Obliczenia własne. 


Jak łatwo zauważyć, statystyka Q pozwala na poziomie istotności 
a = 0,01 uznać za zgodne z danymi empirycznymi modele: (23) i (25). 
Test równoważności widma mocy wskazuje natomiast na zgodność mo- 
delu (25). Obserwacja wyrazów ciągów: {w,T} i (w,F) 17 sugeruje, że ciągi 
te w przypadku modelu (23) znacznie różnią się od siebie. Dużo bliższe 
wydają się być (w;T) i {w,£} dla modelu III. Widać więc, że w omawia- 
nym przypadku statystyka Der? daje bardziej wiarygodne rezultaty. 
Przedstawione powyżej uwagi pozwalają stwierdzić, że optymalny model 
opisujący kształtowanie się cen giełdowych złota w latach 1977—1979 ma 
postać: 


w, = a,—0,19a,_; 
(0,0062) 


ZAKOŃCZENIE 


Zaproponowana w pracy procedura oceny zgodności modelu autoregre- 
syjnego .z danymi empirycznymi dała rezultaty bardziej wiarygodne 
w porównaniu z tymi, które uzyskano przy zastosowaniu testu związane- 
go z ocenami funkcji autokorelacji reszt modelu. Wprawdzie prezento- 
wany przykład empiryczny nie pozwala uogólnić powyższego stwierdze- 
nia, ale sugeruje celowość przeprowadzenia symulacyjnych badań służą- 
cych porównaniu mocy rozważanych testów. Zagadnienie to stanowić bę- 
dzie temat dalszych badań autora. Podkreślenia godny jest również fakt, że 
niniejsza rozprawa potwierdza słuszność prostej procedury identyfika- 
cyjnej modeli typu ARIMA, zaproponowanej przez G. E. P. Boxa 
i G. M. Jenkinsa. Wydaje się, że analogiczne badania można byłoby prze- 


17 Ze względu na długość tych ciągów nie prezentuje się ich w pracy. 
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prowadzić dla cen wolnorynkowych różnorodnych artykułów zarówno 
w warunkach gospodarki kaptialistycznej, jak i socjalistycznej. Uzyskane 
w ten sposób modele byłyby niewątpliwie pomocne w opisie szeregu zja- 
wisk ekonomicznych. 
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PE3IOME 


B cTarbe npeąjaraeTca HOBaA uIpolelypa OLĄEHKM COOTBETCTBMA aBTOpe- 
TpeccMBHbIX Mojelleń WaHHbIM Xe/CTBUTEJILHOCTM. OMIUPMUECKOŃ mpoBepKe No- 
BepraeTcA MEeTOĄ MNeHTMCDMKAHMM aBTOPErTpeCcMBHbIx Monqejeń, IpeXJIOXKeHHBTNA 
Am. 9. II. Boxcom n Ax. M. JIIKeHKUHCOM. DScbcpekT NpoBe1eHHbIX UCCJIEĄOBAHMŃ 
npeĄcTABIEH MOĄEJIbIO, ONKObIBAIOLIEŃ pOpMMPOBAHME eH 3OJIOTa Ha JIOHĄOHCKOŃ 
Oupme B roĄqbi 1977—1979. 


SUMMARY 


The paper presents a suggestion for a new procedure of the checking of the 
goodness of fit of ARIMA models. Box — Jenkins method of the identification of 
autoregressive models is subject to empirical verification. The methods presented 
in the paper are verified on the basis of a model which describes the behaviour of 
Stock — Exchange prices of gold in London in 1977—1979. 


